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Abstract. We compute one of the bulk viscosity coefficients of cold CFL quark 
matter in the temperature regime where the contribution of mesons, quarks and 
gluons to transport phenomena is Boltzmann suppressed. In that regime dissipation 
occurs due to collisions of superfluid phonons, the Goldstone modes associated to 
the spontaneous breaking of baryon symmetry. We first review the hydrodynamics 
of relativistic superfluids, and remind that there are at least three bulk viscosity 
coefficients in these systems. We then compute the bulk viscosity coefficient associated 
to the normal fluid component of the superfluid. In our analysis we use Son's effective 
field theory for the superfluid phonon, amended to include scale breaking effects 
proportional to the square of the strange quark mass m s . We compute the bulk 
viscosity at leading order in the scale breaking parameter, and find that it is dominated 
by collinear splitting and joining processes. The resulting transport coefficient is 
C = 0.011 m^/T, growing at low temperature T until the phonon fluid description 
stops making sense. Our results are relevant to study the rotational properties of a 
compact star formed by CFL quark matter. 
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1. Introduction 

In this paper we present the computation of one of the bulk viscosity coefficient of color- 
flavor locked (CFL) quark matter [1] at low temperature. This work represents a follow 
up of Ref. [2j, where the shear viscosity in the CFL phase was computed in the cold 
regime where the contribution of all the gapped particles (mesons, quarks and gluons) 
is Boltzmann suppressed. 

The density reached in the core of neutrons stars might be so high that all the 
hadrons could be melted into their fundamental constituents. This consideration has 
motivated the studies of QCD at high baryonic density and low temperature [3]. At 
least in the asymptotic high density regime, when the QCD coupling constant is small, 
reliable theoretical predictions for the behavior of quark matter can be formulated. 
Further, it has been known for long time that cold dense quark matter should exhibit 
the phenomenon of color superconductivity. In order to connect theoretical predictions 
with possible astrophysical signatures of quark matter [4], it is important to have a 
precise knowledge of both the equation of state of quark matter and also of all transport 
coefficients, which are very sensitive to the presence of superconductivity in the system. 

It has been established that the viscosities put stringent tests to astrophysical 
models for very rapidly rotating stars, such as for millisecond pulsars. This is based 
on the existence of r(otational)-mode instabilities in all relativistic rotating stars [5], 
which are only suppressed by sufficiently large viscosities. So the viscosities allows 
to discard unrealistic models for millisecond pulsars. There are many different color 
superconducting phases, and their occurrence depends on both the values of the 
baryonic chemical potential and of the different quark masses. At present, there 
are several computations of viscosity coefficients in the different quark matter phases 
[3 [7J [HI [9j [TOl [TlJ [12]. All of them have the motivation of studying the development of 
the r-modes of a compact star made of, partly or entirely, unconfined quark matter. 

Here we will only be concerned about the CFL phase, which is the preferred phase 
in the presence of three light quark flavors. The CFL phase is special in many ways, 
as its long distance physics is very similar to the corresponding one of QCD in vacuo 
[U [13]. Here we want to stress that it is also very peculiar for its hydrodynamics. 
In the CFL case the baryon symmetry is spontaneously broken, and thus CFL quark 
matter becomes a superfluid of the same sort as those found in condensed matter 
systems, such as in Bose-Einstein condensates. Landau developed the hydrodynamical 
description of these non-relativistic superfluids, proposing his famous two-fluid model 
[T4l [15]. The hydrodynamics associated to relativistic superfluids has been much less 
studied, although the two-fluid model has been generalized to the relativistic domain 
[16| \T7\ [181 [19] . We believe that the CFL superfluid may represent one specific and 
beautiful example where the sophisticated relativistic superfluid hydrodynamics could 
be derived from first principles, at least in the asymptotic high density domain. One 
of the peculiarities of these superfluids is that they have more viscosities than a normal 
fluid, as one can define more than one hydrodynamical velocity. 
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In this article we compute the bulk viscosity coefficient associated to the normal 
fluid component of the cold CFL superfluid. Transport coefficients are very sensitive to 
the temperature T of the system. In the regime where T is smaller than all the energy 
gaps of all the quasiparticles (mesons, quarks and gluons), transport coefficients in the 
CFL phase are dominated by the collisions of the superfluid phono ns [201 |2] , and these 
will be the relevant modes in our study. In our analysis we use Son's effective field theory 
for the superfluid phonon [21J. However, Son's theory is scale invariant, and leads to 
a vanishing bulk viscosity, as this transport coefficient measures the dissipation after a 
volume compression or expansion of the system. We then consider scale breaking effects 
due to a non- vanishing value of the strange quark mass. The inclusion of quark masses, 
that represent an explicit chiral symmetry breaking effect in the QCD Lagrangian, 
makes the octet of (pseudo) Goldstone bosons of the CFL phase (the pions, the kaons, 
the eta) massive. Because a quark mass term in the Lagrangian respects the baryon 
symmetry, the superfluid phonon remains massless, although quark mass effects still 
affect its dynamics, as we will see. We compute the bulk viscosity at leading order in 
this scale breaking parameter, and find that it is dominated by collinear splitting and 
joining processes. Surprisingly, the computation shares many technical similarities with 
that of the bulk viscosity in the hot, weakly coupled, phase of QCD at zero chemical 
potential [22], as we will later point out. 

Let us stress that at higher temperatures other quasiparticle modes might be 
relevant for bulk viscosity as well. In Ref. [11] the bulk viscosity due to kaons in the 
CFL phase has been computed. Allowing for flavor changing processes, mediated by 
the electroweak interactions, the bulk viscosity has been computed assuming that the 
relevant processes are those of a neutral kaon decaying into two superfluid phonons, 
and to K ± <-> e + v. However, as found in that reference, for temperatures below 
the energy gap associated to the kaon, 5m, the kaon contribution to bulk viscosity 
is exponentially suppressed ~ e ~ Sm ^ T , as naturally expected. There is an additional 
uncertainty of at what temperatures this suppression is effective, as the value of the 
kaon masses computed in the literature can only be trusted in the asympotic high 
density limit. They are believed to be in the range of the few MeV, or slightly higher 

[231 ES ES ES E3 EE]. 

This paper is structured as follows. In Sec. [2]we recall the hydrodynamical equations 
of a relativistic superfluid. Sec. [3] is devoted to review Son's effective field theory for 
the superfluid phonon. With the same effective field theory, one can easily also derive 
the dispersion law for the phonon in a moving superfluid, introducing the concept of 
acoustic or sonic metric, which also allows us to identify these phonons with the sound 
waves of the superfluid. In Sec. @] we see how scale breaking effects could be included 
in Son's Lagrangian. The explicit computation of the bulk viscosity is given in Sec. HI 
After identifying the leading collisional processes relevant for this transport coefficient, 
we write down the Boltzmann equation for the phonon, and linearize it in the small 
deviations around equilibrium in Subsec. 15.11 The collision term is explicitly written 
down in Subsec. I5.2[ and the numerical results of our computation are displayed in 
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Subsec. 15.31 We conclude with a discussion of our results in Sec. We will use 
throughout natural units h = c = k^ = 1 and the metric conventions (1,-1,— 1,-1). 

2. Hydrodynamics in relativistic superfiuids 

In this Section we present a quick review of the hydrodynamical equations for a 
relativistic superfluid. These represent the natural relativistic generalization of Landau's 
two fluid model of superfluid (non- relativistic) dynamics [T4"l [T5] . There are different 
formulations of the hydrodynamics of a relativistic superfluid [HI [T7J [181 [12], they all 
differ in the choice of the hydrodynamical variables used to describe the fluids. 

The hydrodynamical equations in the superfluid take the form of conservation laws, 
as in a normal fluid. If n p is the particle current (in our case, the baryon current) these 
are 

d p n p = , (1) 
and the energy-momentum conservation law 

d p T pa = . (2) 
In the absence of dissipation, the entropy current is conserved, and thus one further 

has 

d p s p = . (3) 

In the superfluid, the gradient of the phase of the condensate allows to define the 
four vector 

= d p (f . (4) 

The approach of Carter and Khalatnikov [17] defines all the hydrodynamical 
equations based on expressing both the particle current and the energy-momentum 
tensor in terms of \x p (or ip) and s p . In particular, in Ref. [18] it is shown that the 
energy-momentum tensor can be written as 

T pa = Afi p fi a + Bs p s a - Pg pa , (5) 

where P is the generalized pressure of the system. The coefficients A and B can be 
obtained with the knowledge of P. Similarly, n p can be expressed in terms of both /i p 
and s p . 

Son formulated a different description of the hydrodynamics of the relativistic 
superfiuids in Ref. [19] . After a non-trivial mapping of his variables, their equations 
can be converted to the Carter and Khalatnikov form [29] . 

In the zero temperature limit, when the pressure is only a function of the chemical 
potential, the entropy current vanishes. It is only in this case when the the energy- 
momentum tensor takes the form of that of an ideal fluid [18] . If we define the velocity 
vector 

u p = — — (6) 
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such that is properly normalized, u p u p = 1, then in the cold T — > limit one has [TS] 

T pa = (n/j,)u p u a - Pg pa = (e + P) u p u a - Pg pa (7) 

where e is the energy density of the system, and we have used the zero temperature 
relation n/i = e + P. In this case the associated hydrodynamical equations are the same 
as in an ideal fluid. The entropy strictly vanishes, and thus, there is no dissipation. 

At finite temperature, the entropy does not vanish, and dissipational processes are 
responsible for entropy production. Dissipative relativistic superfluid hydrodynamical 
equations have only been derived, to the best of our knowledge, in Ref. p3], although 
there is a vast literature on the subject for non-relativitistic superfluids [HI [15]. In 
Ref. [16] . and after imposing that deviations from the dissipationless particle current and 
energy-momentum are expressed in terms of linear gradients of the basic hydrodynamical 
variables, and that entropy production is positive-definite, it was found that more kinetic 
coefficients that in a normal fluid can be defined. We leave for a future project a much 
more detailed discussion on all the possible transport coefficients that can be defined 
in the CFL superfluid. We simply note here that in a non-relativistic superfluid four 
viscosity coefficients can be defined [TH [T5] , and thus, at a minimum, the same number 
of viscosity coefficients in the relativistic domain exist. 

In this paper we compute the bulk viscosity coefficient associated to the normal 
fluid. For that purpose, and taking into account that the bulk viscosity is a Lorentz 
scalar, we will work in the superfluid rest frame, as this simplifies enormously the 
computation. More specifically, we will compute the dissipative term in the energy- 
momentum tensor that goes as 

r d j = -(5 ij V-V (8) 

where V is the velocity of the normal fluid in the superfluid rest frame. 

3. The cold CFL superfluid and the sonic metric 

In this Section we review the effective field theory of the superfluid phonon constructed 
by Son [21] . Further, we introduce the concept of sonic or acoustic metric, which is rather 
convenient in order to describe the dynamics of the phonon moving in the background 
of the superfluid. It also allows us to give the interpretation of the superfluid phonon 
clS db sound wave. 

In a non-relativistic superfluid, gravity analogues for the description of the low 
energy collective modes or superfluid phonons have been developed (see [3QI ED] and 
references therein). In such an approach, one treats the superfluid as a gravitational 
background, in which the quasiparticles, composing the normal fluid, propagate. We 
will use the same analogy here. 

Son showed that the effective field theory for the only truly Goldstone boson of 
the CFL phase can be constructed from the equation of state (EOS) of normal quark 



Bulk viscosity in the cold CFL superfluid 



6 



matter [21] . If (p is the phase of the condensate, and one defines D p ip = d p (p — (/i, 0), 
then the effective Lagrangian for (p is expressed as 

C eS [D p( p] = P\pL = (D pV D^y/ 2 ] , (9) 

where P is the pressure of the system at zero temperature. 

At asymptotic large densities the EOS of CFL quark matter reads 

m = i-y > (io) 

where /i is the quark chemical potential. At very high /x, when the coupling constant 
is small g(fi) <C 1, the effects of interactions and the effects of Cooper pairing are 
subleading and neglected in Eq. (fTUl) . Also, one assumes that the quark masses give a 
subleading effect, as m q <C /j. From Eq. (flOl) Son obtained the effective Lagrangian 
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There is an interesting interpretation of the equations of motion associated to tp. 
The classical equations of motion derived from the above Lagrangian can be re-written 
as the hydrodynamical conservation law of a current representing baryon number flow, 

d u (n u u ) = , (12) 

where n = j^\ p=po is interpreted as the baryon density [21]. Son defined the superfluid 
velocity u p as being proportional to the gradient of the condensate phase [H], 

u p = , (13) 

where /Jo = (Dp^DPip) 1 / 2 . It only differs from the choice of the last Section by an 
irrelevant constant, although the hydrodynamics is completely analogous. The energy- 
momentum tensor associated to the theory described in Eq. (fTTj) can also be written in 
terms of the velocity defined in Eq. ( TTBl and Noether's energy-density e, 

T pa = (e + P)u p u a - g pa P . (14) 

It is conserved and traceless 

d p T pa = , T p = . (15) 

Eqs. (Tl2]) and (fTBT) are the hydrodynamic equations for the relativistic superfluid 
[21] . They need modifications at finite temperature as explained in the previous 
Section. At low temperatures, the superfluid phonons are thermally excited and conform 
the normal fluid component, which are responsible for the entropy current in the 
system. Other particles may also conform an additional component to the normal 
fluid, but at low temperatures, as discussed in the Introduction, their contribution to 
the hydrodynamics is Boltzmann suppressed. 

Let us also mention that Son's Lagrangian yields the effective field theory of the 
phonons moving in the background of the superfluid. The phonon is a Goldstone 
boson, given also by the phase of the condensate. When the superfluid is at rest, their 
interactions are given in Eq. (fl~TT) . To find the phonon dispersion relation in a moving 
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superfluid we will simply consider the quantum fluctuations around the classical solution 
of Eq. (dU) 

(p(x) = <p(x) + (f)(x) . (16) 
The action associated to Son's Lagrangian 

%>] = Jd'xC eS [d^} (17) 
is then expanded around the classical solution 

The equation of motion of the linearized fluctuation - here the superfluid phonon - 
can be written as that of a boson moving in a non-trivial gravity background 

dJV^gg^d v 4>) = o (19) 



where in this case, we identified 



5 2 CeS 



6(du<p)6(d v <p) 




(20) 

where c s = is the speed of sound in the ultrarelativistic system. Thus, from Son's 
Lagrangian, we have derived the so-called sonic or acoustic metric tensor Q^ v [18, 30J. 
One then finds the phonon dispersion relation in the moving superfluid as solutions of 
the equation 

g^k^k u = o , (21) 

where k^ = (E, k). 

Let us note that in an ideal fluid the fluctuations of the pressure, 5p, obey the 
sound wave equation (see e.g., Eq. (23) of Ref. [32]), which in Fourier space reads, 

. (22) 




u^kf, 



Simply by defining 



2 



^=|^W (23) 

we derive the same dispersion equation (I2T|) . This is naturally so, as the superfluid 
phonons describe the sound waves associated to the superfluid component of the fluid. 
Let us stress here that in the relativistic superfluid there are two (first and second) 
sound speeds, associated to the fact that there are two different fluid components [TH] . 

In the superfluid rest frame, that is, where = (1,0,0,0) the phonon dispersion 
relation simplifies to the form 

E k = c s k , (24) 

as it can be easily checked. This is the frame where we will perform our computations. 
However, we want to stress the fact that this would not be correct in order to compute 
other transport coefficients in the system that involve the superfluid velocity u^. In 
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those cases, the gravity analogues described here become a very efficient tool, which 
will be exploited in future references. Although sonic metrics are usually employed to 
construct simple laboratory analogue models of Einstenian gravity [33], we think the 
method also shows promise in reverse, for the involved features of relativistic superfluids 
might be conceptually understood in the framework of gravity analogues. 

4. Effective field theory for the superfluid phonons including scale breaking 
effects 

In Ref. [2] the effective field theory Lagrangian (TTTT) was used to compute the shear 
viscosity in the phonon fluid of cold CFL quark matter. For the computation of the 
bulk viscosity one has to introduce corrections to that effective field theory. As signaled 
by the vanishing of the trace of the energy- momentum tensor, Eq. ( fTBT) . this is an scale 
invariant theory. Bulk viscosity is a transport coefficient that measures the dissipation 
after a volume compression or expansion, and it vanishes exactly for a scale invariant 
relativistic theory. Thus, in this Section we will look for corrections to Eq. ( TTTT) that 
introduce scale breaking effects. 

The quantum scale anomaly breaks the conformal symmetry in the system 
introducing, through dimensional transmutation, the quantum scale Aq CD . One could 
then compute (^-corrections to Son's effective field theory, arising as ^-corrections to the 
pressure of quark matter. These corrections would introduce new terms in the superfluid 
phonon Lagrangian proportional to the QCD beta function. In the very high /i limit, 
when g(fi) < 1, we expect this to be a rather negligible effect. 

The inclusion of quark mass effects in the system also breaks scale invariance. 
In CFL quark matter, the quark masses represent an explicit chiral symmetry breaking 
effect, which gives masses to the associated (pseudo) Goldstone bosons (the pions, kaons 
and eta). A quark mass term in the QCD Lagrangian respects baryon symmetry, and 
thus it does not make the superfluid phonon massive, however it still affects its dynamics, 
as we show below. 

Since all three light quarks participate in the CFL phase, the largest effect comes 
from the strange quark mass m s , and we ignore que masses of the up and down quarks, 
as m u ,m^ <C m s . Further, we will consider that m 2 < 2A/z [3], which is the threshold 
value under which the CFL phase is stable. We will always work in a leading order 
expansion in m 2 /^ 2 . After imposing the constraints of electrical neutrality and beta 
equilibrium of quark matter, the first correction at order m 2 / '/z 2 to the pressure reads 



To this order, first in the m 2 /fi 2 expansion, and since isospin breaking effects happen 
to be of the same order, one must specify that /i refers to precisely fi n /3 in terms of the 
baryon chemical potential. Let us stress here that as in Son's theory, Eq. (fTUj) . we neglect 
both the effects of interactions and of Cooper pairing, assuming the asymptotic high 
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density, and thus the weak coupling domain. The expansion here is aimed to consider 
the most relevant scale breaking effect that corrects Son's Lagrangian. 

Following the same procedure advocated by Son, with the knowledge of the pressure 
we get the effective field theory for the phonons or sound waves in the system 
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We now rescale the phonon field 

3/i \. 



71 



6/r 



(26) 



(27) 



to normalize the kinetic term in accordance with the LSZ formula. Then the Lagrangian 
for the rescaled field reads 

£ eff = l -(d^) 2 - i(dw) 2 - gMd^d^) + gAd^d^f , (28) 



where, to order m 2 /fi 2 we find 



9s 





m 2 \ 




— ( 


' m 2 

+ v 



TV 



9i 



1 + ^ 



(29) 



(30) 



108/i 4 

Thus, the leading order effect of the strange quark mass in the phonon Lagrangian is 
to modify the velocity of the phonon, that obviously equals to the speed of sound in 
the conformally broken system, and a finite renormalization of the cubic and quartic 
self-couplings. 



5. Transport theory in the phonon fluid 

We now turn to the microscopic description of the bulk viscosity at low temperature. 
The dominant process relevant for the bulk viscosity, shown in Fig. [T], is phonon collinear 
splitting or joining processes 1 <-> 2. Large angle 2^2 scatterings are suppressed, as 
compared to collinear splitting, by powers of l//i 2 as found in Ref.[2]. Because we are 
considering the regime T /i, we can safely neglect those collisions. Further, we don't 
consider small angle 2^2 collisions, which are collinearly enhanced [2]. Considering 
simultaneously the 1 «-> 2 processes and small angle 2 «-> 2 collisions would mean to 
incur in a wrong double-counting. 

It is remarkable that the computation in the asymptotic large density and low 
temperature CFL phase has many points in common with the same computation in the 
very hot, weakly coupled phase of QCD at vanishing chemical potential. In the hot phase 
of QCD, bulk viscosity is dominated by both effective collinear splitting processes 1 «-> 2, 
as well as by 2 <-> 2 collisions [22J. Fortunately, in the CFL phase the computation is 
simpler, as the last processes are certainly suppressed. However, we will find convenient 
to follow the same technical strategy that in Ref. [22J, and tackle different subtle points 
in the computation in the same way as in that reference, which we will closely follow. 
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Figure 1. Under spherically symmetric, radially non- uniform rarefaction 
(compression), of a gas not invariant under dilatations, the pressure diminishes below 
(increases above) its equilibrium value. For the superfluid phonon gas, equilibrium is 
restored by phonon collinear splitting (joining). 



There different technical remark that we would like to point out here. 

Collinear splitting is closed by a convex dispersion relation when one considers 
corrections to the phonon dispersion relation at order k 2 /A 2 [35], where A is the 
superconducting gap. If one considers those corrections, one should then study other 
number changing collisions, such as a 2 <-> 3 scattering for the computation of the bulk 
viscosity. In these 2 «-> 3 processes, one may treat one of the particles as a spectator 
whose only role is to restore energy-momentum conservation in the splitting vertex. We 
will ignore this subtlety here and operate as if the dispersion relation was exactly linear, 
kinematically allowing for collinear splitting and rejoining. We expect that considering 
those processes would simply allow us to get corrections of order T 2 /A 2 to our leading 
result. 



5.1. The Boltzmann equation and bulk viscosity 

In the superfluid rest frame the out of equilibrium phonon distribution function evolves 
according to the Boltzmann equation 

^ + v p -VJ p = -C[f p } , (31) 

where v p = V P E P , and C[f] is the collision integral. We have also introduced the 
short-hand notation f p = f(x, p) that we will use in what follows. 
We will consider small deviations from equilibrium 

fp = f° P q + fl + ■ ■ ■ ■ (32) 
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In the phonon rest frame 

fp q = fp = e p Ep _ i , ( 33 ) 
where j3 = 1/T. 

The Boltzmann equation is linearized in the departures of equilibrium. One has to 
keep in mind that the collision term evaluated with the equilibrium function vanishes by 
detailed balance. As in Ref. [22], we cast the left-hand-side of the Boltzmann equation 
with separated variables as 

X(x)(3f p (l + f p )q(p) (34) 

that parametrizes the advective derivative of / by defining q(p). In order to do so, 
one needs some thermodynamical relations. In the superfluid rest frame, those are 
e + P = Tg, exactly as in the hot quark-gluon plasma at zero chemical potential. The 
explicit calculation of this advective term in [22] yields 



and for a departure of equilibrium that is a uniform compression, 

X(x)=V-V (36) 

as relevant for the bulk viscosity. 

Eq. (1351) reflects the fact that in a relativistic scale invariant theory, where c 2 — |, 
q{p) = and thus the bulk viscosity vanishes. Thus, we allow for scale breaking effects, 
as anticipated in the previous Section. A non-vanishing strange quark mass, m s ^ 0, 
introduces a correction to the speed of sound, as found in Eq. (129|) . In this case we 
obtain 

- if • < 37) 

proportional to our scale breaking parameter m 2 s / 'fj, 2 . 

The next step in the analysis is to project the Boltzmann equation into weak 
(integrated over p) form amenable to variational treatment. To alight notation it is 
convenient to introduce a scalar product 

Multiplying the Boltzmann equation by — T 2 we obtain 

- TfK 1 + /p )zt^*(*) = -T 2 C[f p + fl\ , (39) 

and we define, as in [22], a function for the left-hand-side source S, without the X(x) 
factor (that will cancel left and right in the Boltzmann equation), nor the conformal 
symmetry breaking factor, (m s //j,) 2 , 

S(p) = -Tf° p (l + f° p )^-. (40) 
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The projection of this equation over a complete orthonormal basis of functions ip n is 
then, formally 

2 

Tf) 

X(x)-i^ n \S) = ^ n \C[f p +f p ]) . (41) 

And we further introduce a dimensionless variable (henceforth, the bar notation denotes 
an adimensional function of momenta over temperature) 

S(p) = -T 2 S{p/T) . 

The bulk viscosity is given as the matrix element 

C ^ (S |gZ)-'| S > (42) 
with the differential matrix of the collision operator 

(§) = <*.I^Ww>- (43) 
5.2. The collision term 

Bulk viscosity involves the relaxation of a momentum gradient along the same direction 
of the momentum. Within Son's effective theory, the lowest order effect that achieves 
this is the collinear splitting induced by the cubic term. Because our computation is 
done at leading order in the scale breaking parameter, m 2 //z 2 , and this has already 
taken into account in Eq. ( |35i) . it is enough to keep the scattering matrices as arising in 
the scale invariant theory. Thus, in this Section, we take the value c s = l/v3. Let us 
insist again on the fact that collisions involving two bosons in the initial state and two 
bosons in the final state are suppressed by a further power of /z 2 at the amplitude level, 
and thus we neglect them. 

We will need therefore the amplitudes for a boson of momentum p to split into two 
bosons of momenta p', k' and the amplitude for a boson of momentum p' to split and 
give back in the final state the p boson and a k' boson. From Son's Lagrangian density 
these are found to be, 

M(p; p', k') = (pV • k') + p °(p ■ k') + k'°(p ■ p')) , (44) 

or, employing momentum conservation and the linear dispersion relation p° = c s \p\, 

An 2 
81/i 



\M(p;p',k')\ 2 = —c 2 s \p\ 2 \k>\Hx 2 (\p\ - \k'\) 2 (45) 



where x = p-k' . Similarly, one finds 

4vr 2 



\M(p';p,k')\ 2 = ^c 2 M 2 |A;'| 2 4* 2 (M + \k'\) 2 . (46) 



81/i 
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The form of the collinear splitting collision terms that enter into the Boltzman 
equation can be read off reference 



c L -Vp]=crvp\+crVp) (47) 

1 ' r ]M2n)W A \p-p' -k')(f p (l + f p r)(l + f k ,)- + f p )) 



4E P 



U(27r)V 4 )(p' —p — k') (f p ,(l + f p )(l + fr) - f 9 f v {l + f p ,)) 



J v 

where we introduced the shorthand notation 



2£V(27r) 3 2£y(27r) 3 ' 

Let us reduce Eq. ( 1471) to a more tractable form. The idea is to split the momentum 
conservation delta into an energy part, a longitudinal part (defined along p), and a 
transverse part. The best way to organize the calculation is to introduce a collinear 
splitting function 

i ii n f d 2 k' ± d 2 p' ± \p\ 

\M(p;p',k')\W>( Pl _ -p' ± - k' ± )S(E p - E„ - E v ) , 
so that the first collision term in Eq. (14~TI) becomes 

C l"U] = 2|^/ / d|Wl7(P5l/,W|p| -Pt-^L) (49) 
(f p (l + f p ,)(l + f kl )-f p ,f k ,(l+f p )) 

Now, energy conservation in Eq. (|48p forces the transverse momentum to vanish 

6(c s (\p\ - \p>\ - \k'\)) = , . (50) 



Thus 

7 (p; p' = p - fc', *0 = (^r) W 2 l*?(bl - l^'D 2 • (si) 

Finally one can express the collision term as 

9<7r roo roo 

c'r 2 = ^ / / d\k'\d\ P Up;p',k')5(\p\ - W\ - \*\) (52) 



2\p\ 2 Jo Jo 
{f P {l + f P r){l + fv)-f P rfk>(l + f P )) , 

with 

7W=P-^V) = |^jblWl J . (53) 

We can analogously reduce the second term of Eq. (jl7j) . that amounts essentially to 
the exchange p —> p' with respect to the C^ 2 piece, as the splitting function is totally 



/ P 1 = ^?/ P °(l + f P )^- (54) 
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symmetric in its three arguments. Thus, we won't write here the analogous equation 
for C^ 2 . 

We linearize the collision term in the first order in the deviation from equilibrium, 
which we parametrize as 

The last line in Eq. (I52"|) containing the distribution functions becomes 

F(p; P ',k>) ee X jf )Xp, ~y~ Xp / p °(l + ff)(l + ^) • (55) 

We then expand the deviation from equilibrium in terms of the function basis 
xip)/T = J2 n Xn^nip/T). Then the projected, linearized collision operator turns out to 
be 

(iM^M = (56) 

where we define F n as the value of F evaluated at tp n , instead of at x/T- 

To obtain the parametric dependence of the bulk viscosity we define the 
dimensionless quantities 

7 = ^7 (57) 

in terms of which, and extracting the factor X(x)/T, the right hand side of Eq. ( |56l ) 
turns into 

X{x)T 5 5C_, 

where the energy-dimension is explicit, since C mn is now a function of the ratios of 
momenta over temperature alone. One can then check the dimension of the projected 
Boltzmann equation 

^ :S{P) = J^W) (59) 



-(^mlTTl^n) (5* 



(2vrT) 3 u> /i 4 5/ 

that matches the defining Eq. (j4"0l above. 

Finally, the parametric dependence of the viscosity, following from Eq. (|42j) . is given 

as 

C = ^l(§)~V>- (60) 
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5.3. Numerical evaluation 

Once the parametric dependence of the bulk viscosity is known, all that remains is 
to evaluate a numerical factor. A subtle point comes in choosing an appropriate trial 
function family ipi, so that all integrals converge appropriately in both their infrared 
(IR) and ultraviolet (UV) domains. For example, the natural family of orthonormal 
functions in the interval (0, oo), the Laguerre functions, would yield 

^ m (p) = ^ L m (p) (61) 

with conventional L m Laguerre polynomials. The 1/p 2 from two such functions, together 
with a 1/p from the Bose-Einstein factor (e p — 1) _1 — > p^ 1 and our splitting function 
7 oc p 2 , make the integral in Eq. (1561) infrared divergent. Fortunately, it is not necessary 
to use an orthonormal function family 0, as long as one is not interested in the function 
X itself, but only on its projection to obtain the transport coefficient. This can be shown 
with minimum linear algebra by changing basis from an orthonormal to an arbitrary ipi 
family function. We therefore choose 

(p/T) m 

^(P/T) = (1 ■ ( 62 ) 

One should also notice that there is a zero-mode of the collision integral, visible in 
the last line of Eq. (1551) . that vanishes when ip — p for p = p' + k'. The collision operator 
is not invertible, but as observed in [22], this does not suppose a problem as one can 
add to C an arbitrary constant A times the projector over this zero mode 

C ^C + \\E p p 2 )(E pP 2 \ 

since its projection over the source in Eq. (!60|) vanishes. Note a factor of + /°) can 
be multiplied to the vector \E p p 2 ) without numerically affecting the result, as in Ref. 
[22]. 

In table [1] we show the fast convergence with the number of functions employed 
(size of the linear system). The integration is performed with a Gaussian grid, there is 
no sensitivity to UV or IR cutoffs, as shown in Fig. [2], nor to the parameter A that fixes 
the zero-mode subtraction. In the table however we give two sets of numbers showing 
that the result is essentially equivalent should the factor / (1 + fo) be omitted in the 
zero-mode subtraction. 

6. Discussion 

We have established that the bulk viscosity coefficient associated to the normal fluid 
component of a cold CFL superfluid is given by 

4 

777, 

Ccfl = O.Oll^f (63) 
| We thank Guy Moore for this observation 
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Table 1. Convergence of the pre-coefficient of £ with the number of functions. 
First column: number of functions in the family. Second column: result with zero- 
mode subtraction based on the function E p p 2 . Third column: result with zero-mode 
subtraction based on the function E p p 2 f p (l + f p ). This pre-coefficient can be thought 



m s 


= 1,T=1) 




m 


Cl 


(2 


1 


0.00978 


0.0093 


3 


0.0109 


0.01057 


5 


0.0110 


0.01067 


7 


0.0110 


0.01070 




Figure 2. Integrand for the last integral in Eq. (|56|) yielding the dimensionless £, 
showing IR and UV integrability. 

at first order in the conformal-breaking parameter m^/fi 2 . This is a remarkable result 
for several reasons. First, it is much smaller than the shear viscosity already reported 
in Ref . [2] . The reason for being so is that the two coefficients are governed by different 
sort of processes (collinear splitting for bulk viscosity, large-angle collisions for shear 
viscosity), which occur at different rates. Second, at the order we computed it is seen to 
be independent of the chemical potential and the superconducting gap A. Third, due 
to the dynamics being dominated by the superfluid phonon, which is a Goldstone boson 
which remains always massless, it is not exponentially suppressed as might be thought 
of based on a calculation involving gapped degrees of freedom, as, for example, those 
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due to quarks [8] or to kaons [TT|. Thus, it is the leading contribution at very small 
temperatures T <C m, A, where m is the energy gap associated to the lightest massive 
mode of the CFL phase. This is the temperature regime that we call "cold" . 

The result in Eq. ( |63|) could have been anticipated from the mean free path for 
small angle collisions discussed in our previous work [2], namely 

Asmall OC — . (64) 

A quick estimate of the parametric dependence of the bulk viscosity would be 

4 

771 

C oc A smaU x n x (p) x C 2 oc y (65) 

in terms of the phonon number density ri oc T 3 and average momentum (p) oc T, and 
of the conformal breaking parameter C = (~ — c^j ~ m^//i 2 . This immediately yields 



Eq. (16311 up to the numerical factor. Note that the shear viscosity however has a different 
parametric behavior with the temperature, as small angle collisions are very inefficient 
for transferring transverse momentum. This is seen in the calculation in [2j by the 
appearance of near-zero modes that appear in the computation of the 2 —>■ 2 collision 
operator and make collinear splitting irrelevant there [§|. 

While our bulk viscosity result diverges in the limit T — > 0, it should be kept 
in mind [2] that when the temperature diminishes the mean free path of the phonon 
becomes large, and at some point the hydrodynamical description of the phonon fluid is 
meaningless (one rather has free streaming of phonons) . Then only the perfect superfluid 
with no dissipation remains. For astrophysical applications, and assuming that the 
radius of the compact star is of the order of R ~ 10 km, this happens at T ~ 0.06 MeV 

0- 

Although we have computed the bulk viscosity in a fluid at rest, it is possible to 
extract from our results the frequency dependent bulk viscosity needed for astrophysical 
applications. Based on our results, a quick estimate of a so defined frequency- dependent 
bulk viscosity has recently appeared [38]. The results are very interesting and suggest 
that the rate of equilibration of bulk distortions in a hypothetical quark star, at 
physically relevant frequencies, receive contributions not only from weak-equilibration 
processes but also from the phonon-splitting processes that we have studied here [38J. 
Although quark-gluon equilibration times are short, the phono n-system is described by 
the weakly coupled effective Lagrangian of Son, and the bulk viscosity is proportional to 
the (small) scale-violating parameter. This corrects one's first intuition about neglecting 
strong interaction phenomenology altogether in the belief that strong interactions should 
permanently be in equilibrium. 

There is also another subtle point. In the existing literature where one needs a bulk 
viscosity coefficient to perform the analysis of the fate of the r-modes, the computation 

§ The exact zero mode appearing in Eq. (|55p however, is not a separation from equilibrium as it 
maintains the detailed balance relation C[/ + zero] = and can be subtracted. This explains the 
failure of relations one could have guessed such as C oc m*rj/n 4 , that seem to fail analogously in <f> 
theory [37] but hold in the quark-gluon plasma |22j . 
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is performed by analyzing the energy dissipated after a compression or rarefaction over 
one period, r = 2tt/uj, where u is the frequency of the fluctuation. Thus 

^H^CV-V)', (66) 

where V is the hydrodynamical velocity. The resulting bulk viscosity coefficient is then 
given also as a function of u. While this relation is valid for a normal fluid, it should 
be generalized for a relativistic superfluid. In a superfluid there are at least three 
bulk viscosities, and unless the remaining coefficients vanish for CFL quark matter, 
they should contribute to dissipation in an oscillatory compression or rarefaction of the 
system, and should contribute to the right hand side of Eq. ( 1661) . 

It is thus urgent a computation of the remaining viscosities of the CFL superfluid, as 
well as a careful study of its low temperature hydrodynamics. Both are required for the 
study of the r-modes of a hypothetical compact star made of CFL quark matter. Let us 
point out that the relevance of the existence of several viscosities in a superfluid neutron 
star has only been emphasized in very recent publications [3HJ HOI E] • In particular, 
only in Ref. [H], all the bulk viscosity coefficients in a neutron superfluid have been 
computed. 
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